Thermoelectric properties of junctions between metal and strongly correlated 

semiconductor 
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We propose a junction of metal and rare-earth compound semiconductor as the basis for a possible 
efficient low-temperature thermoelectric device. If an overlayer of rare earth atoms differing from 
the bulk is placed at the interface, very high values of the figure of merit ZT can be reached at 
low temperature. This is due to sharp variation of the transmission coefficient of carriers across 
the junction at a narrow energy range, which is intrinsically linked to the localized character of the 
overlayer /-orbital. 



New thermoelectric coolers^and power generators are 
attracting increasing interest E Their quality is governed 
by the dimensionless figure of merit ZT , a product of 
transport coefficients to be defined below times the ab- 
solute temperature. At present, the best ZT obtainable 
is ~ 1 at room temperature, while a value about 3 '--^ 4 
is considered a real breakthrough. 

There are two theoretical ideas for increasing ZT . One 
is to utilize the sharp energy features in bulk jmateri- 
als such as strongly correlated semiconductors. Ha rAn- 
other is to exploit the interface energy structuresQ in 
metal/supercoriductor junctions, which have been tested 
exp er iment ally.El 

In this Letter, we combine these two ideas in explor- 
ing the thermopower behavior of a junction between a 
metal and a particular class of mixed- valent semiconduc- 
tors with high dielectric constants, such as SmBg and 
Sm2Se3, called the electronic ferroelectrics (FE).Q Vari- 
ous cases of the junction are considered, including clean 
and "dirty" interfaces. When a suitable rare-earth impu- 
rity layer forms the interface, we find very high values of 
ZT at low temperatures. There is great interest in low 
temperature thermoelectrics for special applications. 

Mixed-valent semiconductors are rare-earth com- 
pounds, usually cubic, whose relevant electronic proper- 
ties may be modeled by a /-flat band and a broad conduc- 
tion band, with two electrons per unit cell.l3 Strong cor- 
relation between d-electrons and /-holes can renormalize 
the bands and create a temperature dependent small gap 
[see Fig. ^b)] : we consider this case for the semiconduc- 
tor on one side of the junction, and describe it with the 
self-consistent paean field solution (MP) of the Falicov- 
Kimball modeLa The ground state of the insulating phase 
is found to be a coherent condensate of d-electron and /- 
hole pairs, giving a net built-in macroscopic polarization 
which breaks the crystal uiversion symmetry and makes 
the material ferroelectricH 

While previous calculations have dealt with 
metal/ordinary-semiconductor junctions neglecting the 
dependence of the transmission coefficient on the carrier 



energy,H'll3 for the metal/FE junction this feature turns 
out to be crucial, requiring a careful treatment. 
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FIG. 1. (a) Pictorial representation of the junction along 
the z axis perpendicular to the interface, (b) "Bare" band 
structure (A = 0) of materials on both sides of the junction 
(solid lines), and renormalized bands (dashed lines) at T = 0. 
Parameters are t' = 5i and A(T = 0) = 0.49i. 

To this aim, we model the motion along the z direction 
perpendicular to the interface with the one dimensional 
spinless hamiltonian Ti, given by the sum of three terms 
describing the left bulk metal, the interface layer, and 
the right bulk FE, respectively [see Fig. ^a)]: 



Ti = Tin 
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'Wmetai IS a tight-binding hamiltonian for lattice sites at 
z = aj , j < (a = lattice constant), with energy e^ and 
hopping parameter t. T^fe is the bulk Falicov-Kimball 
hamiltonian: 

Hfe = £d XI 4^J + ^/ X! /J+1/2/J+1/2 - < 51 4^J+i 
]>o j>a j>o 

i>o j>o 

+UafY, 4^4+1/2^+1/2 + Uai J2 d]djfU/2f^-y^- (2) 

Here, while electrons can tunnel between d-sites with 
hopping coefficient t <^ t', /-electrons are completely 
localized [/j+1/2 (dj) destroys a /- ((i-)electron at z = 
(j -I- l/2)a (ja) with energy if {id)], but interact with 
d-electrons via the Coulomb repulsion energy Udf- More- 
over, the d-f hybridization Vdf is odd under inversion. 
On this side of the junction we assume one electron per 
unit cell. The hamiltonian term 



'^interface — Sdododo + £/l/2/i/2/l 



/2JI/2 



(3) 



accounts for the possibility of impurity atomic layers at 
the interface, namely one at z = made of d-like sites 
with energy ido and another one at z — a/2 made of 
/-sites with localized orbitals of energy £/i/2- In our 
model, the electrostatic potential V across the junction 
(by applying an external bias) is constant in the two bulk 
regions (the small gap FE semiconductor has carriers at 
finite temperature) and has a sharp step at the inter- 
face. We adopt the following strategy to solve the motion 
across the junction: (i) We apply the MP approximation 
to the hamiltonian H -^ T-C^^ , introducing the built-in 
coherence of the FE bulk A — Udf<. rfI/j±i/2 > (the vac- 
uum is the state with no /-holes) p treating renormalized 
site energies {id ~* £d, £/ ~> £/) as material parameters, 
and choosing e^ = e^ = £/ = 0, namely the middle of 
the d-band on both sides of the junction and the flat /- 
band are aligned [see Fig. |l|(b)]. (ii) We calculate-the site 
coefficients for the Bogoliubov-Valatin operatoiilil 



Ike = 



/N. 



J2 {Mj) d] + Vk{j + 1/2) fl,/,} (4) 



{Ns = number of d-sites) which creates elementary elec- 
tronic excitations of energy w(fc) > (referred to the 
chemical potential fi) if applied to the ground state, i.e. 



ih-jke = [7fce,H^^- fJ-N] = Uj{k)^ke, 



(5) 



where N is th&|number operator, and fc is a suitable quan- 
tum number JIj Similar expressions occur for holes. From 
Eqs. (|j) and (g) one derives equations for the amplitudes 
of the quasi-particle excitation (wfe(j), Vk{j + 1/2)). 

Solutions are travelling waves partly refiected and 
partly transmitted at the interface. In particular, the 
incoming and refiected wave (^il,^2l) is a compatible 



solution for j < plus Uk{0) = *il (0) if *il(J) = 
exp(igaj) - Rkexp{~iqaj), *2l(J + 1/2) = 0, uj{q) = 
eV — 2i' cos {qa) — fi, namely the tight-binding solution 
of the left side metal. Similarly, the transmitted wave 



*ifl(j) 
*2fl(j + 1/2) 



n 



Uk 



Vk& 



ika/2 



^ikaj 



(6) 



is a solution of the bulk FE for j > only if uj{k) = 
£,k + Ek - Mi with S,k = efc/2, Sk = ~2t cos {ka), Ek = 
[^l + |Afc|2 + \Vk\^Y'^, where A^ = 2Acos(A:a/2) and 
Vk — 2iVdf sin {ka/2) have even and odd parity, respec- 
tively, and 2|ufe|2 = 1 + a/-Efe, l"fc|' + l^fcl" = 1- Fig- 
ure |l|(b) shows the quasi-particle band structure on both 
sides of the junction when Vdf — > 0. Note that the FE 
gap is indirect and that the bottom of conduction band 
is much flatter than the top of valence band. 
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FIG. 2. Total (electron plus hole) transmission coefficient 
T(e) +T(— e) vs energy, for different values of impurity levels 
Edo (e/i/2)- Left panel: d-type impurity. Right panel: /-like 
impurity. 

We now have to match left and right bulk solutions at 
the interface. Because quantities 1)^(1/2), T^, and Rk are 
still unknown, the finite-difference equations of motion 
coming from Eqs. (^j and (^ not yet employed are salved, 
with the elastic scattering condition ^(g) — Ci;(fc).ll3 The 
reflection coefficient Tl{ui), defined as the ratio of re- 
fiected to incident fiux, is simply \Rk\ (g, fc > 0), and 
the transmission coefficient is T{lu) = 1 — TZ{u!). 

Electric and heat density currents are calculated to 
first order in the differences of V and T on the two sides of 
the junction. From these the conductance G = 2e'^Lo/h, 
the thermal conductance Gt = '^{L2 — L\/ Lf)) / {Th) , and 
the thepHopower Q = Li/{eTLo) of the interface are 
derived JU where 



L. 



de{T{e)+7-{-e)}e'' 



df{e) 
de 



(7) 



f{x) is the Fermi distribution function, and T{x) ['T{x)] 
is the electron (hole) transmission coefficient, non-zero 
only for positive excitation energies within the band 
range. The interface figure of merit is ZT = Q'^GT/Gt- 
Figures shows the total transmission coefficient T{e) + 
T(— e) vs energy for different values of impurity levels Edo 
(left panel) and £/i/2 (right panel). We consider a sys- 
tem at T = with the same parameters of Fig. |l|(b), 
i.e. one in which the metal band is much broader than 
the semiconductor one and the FE band gap is 0.4t, ap- 
proximately one tenth of the total bandwidth {Vdf — > 0). 
One sees that the clean interface {sdo — £/i/2 — 0) al- 
ready presents a strong electron-hole asymmetry. Hence, 
from Eq. (M) high values of Q follow. The effect of an 
impurity d-layer is just to uniformly depress T -|- T: the 
greater the impurity level energy, the lower the trans- 
mission (left panel). This trend, whose global effect is 
found to enhance ZT, has intrinsic physical limitations 
because Sdo must be of the same order of magnitude as 
bulk energies in realistic systems. The behavior of an 
/-impurity is drastically different. We see that (right 
panel), as we set e/1/2 to negative values, T goes to zero 
in the neighborhood of the same values, as if the hole 
were completely backscattered from the interface when 
resonating with the impurity atom. A similar behavior 
occurs for e/1/2 > 0- The overall effect is so strong even 
to change the dominant (electron or hole) character of 
transport and, hence, the sign of Q. 
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FIG. 3. Two-dimensional contour plot of ZT/{ZT + 1) 
in the T'-e/i/2 plane. The gray scale goes from (white) to 
1 (black), and the contour lines are incremented by 0.1. Pa- 
rameters as in Fig. hi. 

The effect of the /-impurity layer on ZT is shown in 
Fig. H We represent the quantity ZT/{ZT+ 1) (due to 
the huge range of ZT) as a contour plot in the T-efi/2 
plane. While ZT has no upper bound, ZT/{ZT+1) < 1. 
For a fixed value of e/1/2, ^T goes to zero as T ap- 
proaches the critical temperature Tc [KbTc = 0.19i here) 
at which the gap vanishes and the bulk FE turns into a 



metal [T(e) = T(-e)]. In the opposite limit T -> 0, ZT 
can reach very high values (the darkest region enclosed 
by the inner contour line corresponds to ZT > 90), de- 
pending on e/1/2- At r > Tc ZT goes to zero, as e/1/2 
varies, when Q changes sign. 

From Fig. 0, there is an optimal value of e/1/2 of th^ 
interafce layer at each temperature. The corresponding 
maxmimum figure of merit ZT increases with lower tem- 
perature. For example, with a gap of 100 meV at T = 0, 
the best value of ZT will be ~ 1 at T =300 K, but aheady 
- 6 at 150 K and - 100 at 40 K. 

In summary, we have made a qualitative theoretical 
study of the possibility that a junction of metal and 
FE (as opposed to bulk materials) produces high ther- 
mopower. The figure of merit ZT attained is very high, 
especially at low T. In these regimes, bulk thermal con- 
ductivity would be dominated by phonons which would 
reduce ZT. However, if we choose two materials with 
large thermal impedance mismatch, the phonon scatter- 
ings at the interface would decrease the junction thermal 
conductivity.Ej Thus, phonon conductivity would not vi- 
tiate the high ZT found. In a longer paper, we will con- 
sider the role of phonons in greater detail and include 
a more comprehensive study of junctions with different 
classes of semiconducting materials. 
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